How the HBT-Puzzle at RHIC might dissipate 
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I compute the first-order corrections to two-particle Bose-Einstein correlations due to deviations 
from equilibrium. Based on that result 1 argue that for nearly perfect fluids, the HBT radii "freeze 
out" much later than the single-inclusive distribution. Therefore, to prevent a big longitudinal ho- 
mogeneity length, the QGP presumably should hadronize not into a nearly perfect fluid of hadrons 
but rather into a very dissipative one. This could be achieved by a hadronization phase transition 
away from equilibrium. 



Pion interferometry has become a powerful tool for studying the size and duration of particle production in high- 
energy collisions, from e + e~ over pp or pp to heavy ions like Au + Au or Pb + Pb For the case of nuclear 
collisions, the interest mainly focuses on the possible transient formation of a deconfined state of matter. This could 
affect the size of the region from where the pions are emitted as well as the time for particle production. In particular, 
it was hoped that a strong first-order QCD phase transition would lead to long lifetimes of the particle source |2|,|3| . 

Data [Q from BNL-RHIC, and also from the CERN-SPS, represents a big "puzzle" || in that the duration of 
particle emission, as well as the longitudinal homogeneity length, appear rather short. In fact, within the framework 
of relativistic perfect fluid dynamics, they seem incompatible with even a mere cross-over H||. Moreover, studies of 
hadronic rescattering following hadronization of a QGP fluid indicated that one should always observe a long emission 
duration for hadrons if they are formed at the confinement temperature T c ~ 180 ± 20 MeV in equilibrium [Q. This 
' is because the lifetime of the particle source is enhanced by dissipative effects, 
psj , In these notes, I compute the corrections to the single-inclusive pion distribution and to the two-pion correlation 
function due to deviations from perfect equilibrium. I use a simple schematic model which allows for an analytical 
| treatment: one-dimensional (longitudinal) boostinvariant expansion; of course, I must focu s on the longitudinal 



"radius" Rn then. The analytical expression for the longitudinal homogeneity length ( 3.12 ) with first-order non- 
equilibrium corrections is a new result (to my knowledge), generalizing the expression of Makhlin-Sinyukov Q which 
applies to perfect local equilibrium. Based on those results, I discuss why HBT correlation functions are much more 
sensitive to the late, dissipative stages of the evolution than the single-inclusive distributions (in agreement with 
quantitative kinetic theory solutions 0): the first-order theory applied here predicts that the correlation function at 
small relative momentum of the pions can not in fact "freeze out" . This supports the expectation that one should 
observe a long duration of emission if the post-QGP evolution starts with a hot (T ~ T c ) and dense hadron fluid close 
to equilibrium. Hence, the way to go might be: start low. A phase transition away from equilibrium can produce cool 
hadrons with T < T c , and induce large gradients of the temperature and velocity fields. See the last section. 

I do not review here the extensive literature on HBT interferometry for nuclear collisions. Recent studies devoted 
specifically to effects from a confinement transition in Au + Au at RHIC include computations using 

1. ideal hydrodynamics P,pp|-pT||, 

2. ideal hydrodynamics up to hadronization followed by microscopic hadron kinetic theory M, 

3. parton kinetic theory with Jl2] | or without |l5| ] hadronic rescattering, 



4. a QGP evaporation model without interactions among hadrons 14 
Those papers can also be used to trace earlier literature. 



I. FIRST-ORDER DEVIATION FROM EQUILIBRIUM 

A collection of many particles in a large, static volume approaches a steady distribution in momentum space at times 
much larger than the microscopic relaxation time for collisions r c . This is the equilibrium distribution / oq (_E), which is 
either a Bose or a Fermi distribution, depending on the nature of the particles. The temperature T characterizes how 
rapidly the probability to find a particle with energy between E and E + dE falls off as E increases. In equilibrium, 
T is constant throughout the entire volume. 
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More generally, one can think of situations where T varies in space. That is, on scales much larger than the 
interparticle distance but much smaller than the total volume, the distribution of particle energies is still described 
approximately by the equilibrium distribution. If this is not the same distribution everywhere in space, one expects 
hydrodynamic flow to develop between those regions. 

In the absence of a uniform temperature and flow velocity field (global equilibrium), there must also be local 
corrections to the equilibrium distribution. They can be obtained from the Boltzmann equation |l5| ] 



p-0/(x",p")=C[/] 



(1.1) 



The collision kernel on the right-hand-side is a functional of the distribution function f(x tJl 1 p tl ). It vanishes for 
f{x^,p^) = f Qq (E;T(x^)). However, the left-hand-side does not vanish if f(x^,p^) = f eq (E; T(x tl )) with a non- 
uniform temperature field T(x) because of the presence of gradients. Therefore, f cq {E\ T(x M )) is not a solution of the 
Boltzmann equation (1.1). Expanding the collision kernel to first order in deviations from equilibrium, one can define 
the relaxation time for collisions via 



E 

T, 



sc 

Sf 



(1.2) 



Note that r c depends on both momentum and space-time, since 5C/Sf is evaluated at f eq (x^,p^). For simplicity, 
in what follows I neglect the dependence on energy-momentum (more precisely, that on rapidity; p t dependence is 
allowed). 

The local distribution function is then given by 



f(x*,i?) = f eq (E(x^;T(x^)+6f(E(x^);T(x^) ee f^E(x*);T(x»)) l+8f(E(x*);T(x*)) 



(1.3) 



where E(x fl ) = p ■ u{x^) is the energy of the particles measured in the local rest frame. The Boltzmann equation then 
reads 



P ■ d (f cq + Sf) = C [f cq + Sf]~C [/ cq ] + Sf 



(1.4) 



The first term on the right-hand-side vanishes by definition of / cq , while the second term on the left-hand-side can be 
dropped if considering only corrections to first order in gradients Thus, using (1.2), 



p ■ u 



-P ■ df c , 



(1.5) 



In what follows, I focus on the limit m 2 = p\ + m 2 3> T 2 , i.e. heavy and/or high-momentum particles, such that 
the quantum mechanical distribution functions can be approximated by the classical Boltzmann distribution. I also 
assume that no conserved charges such as baryon number are present; in this case viscous corrections to the energy- 
momentum tensor obtained from 8f vanish in the frame where w M = (1,0), which therefore corresponds to the 
Landau-Lifshitz definition of the local rest frame. Then, 



p ■ w(x M ) 



p ■ d 



p ■ u(x^) 



(1.6) 



Below, I will apply (1.6) to compute corrections to the two-particle correlation function. The Israel-Stewart second 
order theory of imperfect fluids in principle represents a more satisfactory, yet much more involved approach |16|. 



II. SINGLE-PARTICLE DISTRIBUTION 



The single- inclusive distribution of particles in a fluid, measured on some space-time hypersurface cr M is given by |T 

dN 



N(p) 



(fptdy 



da-pfip-u) 



(2.1) 
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Consider a fluid with infinite extent in the transverse directions, and with space-like hypersurfaces of homogeneity. 
In particular, for longitudinal scaling flow|^] and vanishing transverse flow velocity those are surfaces of fixed proper 
time t = \ft 2 — x 2 , which is invariant under Lorentz boosts. Then, 

da ■ p — d 2 r dr\ Tm t cosh(y — rj) . (2-2) 

Thus, 

N(p) = / d 2 r drj rm t cosh(y - 77) {/ oq (p ■ u) + Sf (p ■ u)} 

'dWm,co S M»-,ri/„ q (p-u){l + S7(»,,0}. (2.3) 

At large mt/T the Bose or Fermi distributions approach a Boltzmann distribution. Here T denotes the temperature 
on the hypersurface specified by r; for simplicity, let it be independent of the rapidity of the flow. Moreover, in this 
limit the single-inclusive distribution can be evaluated by a saddle-point integration, 



/ oq (p • u) ~ e~^ T = e -"Hcos H y-v)/T ^ _ v ) e -™ t /T 

V m t 



(2.4) 



In other words, evaluate the integral over the flow rapidity 77 by expanding (in the exponential) p ■ u — m t cos h 9 ~ 
m t (l + 2 /2) up to first nontrivial order in 8 = y — rj, and set y = rj throughout the rest of the integrand in ( |2.3| ); this 
reduces the integral over rj to a Gaussian integral. 

The saddle-point approximation also greatly simplifies Sf: 



P ' ^T^a) = j^pjP ' dp ' u ( a "> + P ■ u ( a )P ■ d -f = ~ ^ sinh 2 6> + 77^ cosh ( 



cosh 8d T + — sinh 0cL 

T 



1 



T(a) 



and so 



- rn 2 d T ^ , (2.5) 

Sf(p ■ u) -> T c m t d T ^; . (2.6) 

The integral over d 2 r is trivial, since we assumed that neither the temperature T(a) nor the flow velocity u^^a) 
depend on r; thus, J d 2 r equals the transverse area S t = ttR 2 . We obtain 

^^.^^(l-^^l . (2.7) 

L T 1 OlOgT J 

To be more explicit, we may for example assume that the temperature prior to freeze-out decreases like a power-law 
in 1/t: 

(2.8) 

For ideal isentropic expansion the entropy in a comoving volume element is conserved, which for ultrarelativistic 
particles (m <C T) leads to 7 = 1/3. On the other hand, isoergic expansion is perhaps more realistic near freeze-out; 
in that case the energy per comoving volume element is conserved, which gives 7 = 1/4. Then, 

N{p) = S t T^T^ t e- mt l T [1+^^] . (2.9) 



The leading term in eq. (2.9), corresponding to t c — * 0, is the result for an equilibrium distribution of the particles. 
The second term in the curly brackets is due to corrections from local equilibrium, proportional to the microscopic 



1 This simply means that the flow rapidity r)f = Artanhti^ equals the space-time rapidity 77 = Artanh(x/t) everywhere in the 
forward light-cone, where x denotes the longitudinal direction. 
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relaxation time r c relative to the expansion time r. It gives us a hint about the meaning of "freeze-out" . The hadrons 
measured in the detector are cold T = hadrons; they do not stop interacting suddenly at some temperature^]. 
However, as T decreases in the course of the evolution, corrections to / oq become more and more important. The 
more of them I include in my computation, the lower I can go in T. With the full collision kernel I could go down 
smoothly to T = pq| without introducing an abrupt "decoupling" . The important point to realize though is that as 
corrections to equilibrium grow, they tend to "stabilize" the single-inclusive distribution: from eq. ( |2.9| ) one sees that 
the ideal N(p) at any fixed m t decreases with T, but the correction term acts to slow down that evolution. Eventually, 
at some time r or temperature T, the correction terms "freeze" N(p), i.e. it stops evolving any further despite the 
ongoing approach of T to zero. I can view that as the "freeze-out" temperature. 

What this discussion points at, as well, is that "freeze-out" is not universal but depends on the observable under 
consideration. This is because corrections from deviation from perfect equilibrium are not the same for all observables. 
We shall return to that below. 

The correction to first order in the relaxation time to equilibrium, r c , does not make for a quantitative description 
of the freeze-out process. Nevertheless, to get a rough idea, one might ask when dN(p)/dr = 0. This leads to 
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(2.10) 



Here, To is the time when corrections to perfect equilibrium set in. This could be as early as hadronization, or any 
time after that. To denotes the temperature at that time (mt/To must be large enough so that r c > 0). 

At To we may, for example, be dealing with a perfect fluid, i.e. t c (jo) — 0. As long as t c /t remains negligibly small, 
and t/to is not too big, the solution of ( [2.1C ) is approximately given by 



T - T 



TqTq 1 - 7 /2 
mt 7(1 - 7) 



1-7 
- I I 

TO 



(2.11) 



Thus, t c should grow approximately linearly with the expansion time. As r c becomes larger and comparable to r, its 
growth must speed up in order that N(p) remains "frozen" (also, if it is non-negligible already at hadronization, i.e. 
if the hadron fluid is rather viscous from the start) . To see this, keep only terms proportional to r c on the right-hand 
side of ( 2.10 ), which leads to 



Tc(T) 

t c (t ) 



-7/2 



mt 

ex Pl - 



-1-1 

TO 




— . /__ p m t /T-m t /T 



(2.12) 



Clearly, when the rapid growth kicks in, this signals the break-down of the first-order theory and that corrections of 
higher order in t c also become important. 

So, to summarize, the basic picture is as follows. A hadron "fluid" is produced by the decay of the deconfincd 
state. In theory it may evolve as an ideal fluid for some time, if the relaxation time to equilibrium is extremely short. 
Eventually, as the fluid becomes more dilute, the relaxation time must grow: at first rather slowly, then turning to 
a rapid non-linear growth. Alternatively, the relaxation time may be non-negligible already at hadroniz ation , and 
continue to grow rapidly as the fluid expands. In cither case, as soon as the growth rate of r c satisfies eq. (2.1C), the 
single-particle distribution doesn't change any more, it "freezes" . Smoothly, the fluid cools to T = On. 



III. TWO-PARTICLE CORRELATION FUNCTION 



The two-particle correlation function for identical particles is given by 



2 For example, for a chiral symmetry which is broken spontaneously (by choice of a vacuum) but not explicitly, at low T 
corrections to the pressure of non-interacting pions are determined by chiral perturbation theory: ~ (T//^) 4 pid. 

3 At some point number-changing reactions must also freeze. One then has to introduce a continuity equation for the number- 
current, which will lead to a chemical potential mM. 
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R( Pl ,p 2 ) = N{ Pl )N{ P2 ) [C 2 {p u p 2 ) - 1] 

= Re J da i ■ Kda 2 ■ K exp i(o\ — a 2 ) ■ (pi — p 2 ) 

x {/ oq («i -K)+Sf K • A 1 )} {/ oq (u 2 -K) + Sf (u 2 ■ K)} 



(3.1) 



with = (p^ +p^)/2. I consider particle pairs emitted at the same azimuthal angle and with the same transverse 
momentum, p ti \ — pt, 2 - This is no severe restriction as for 1+ld longitudinal expansion only differences in the 
longitudinal momenta of the emitted particles matter. Introducing the variables a = y 2 — yi for the relative rapidity 
of the emitted particles and 6\ = y\ — r\\ , 8 2 = y 2 — rj 2 for their rapidities relative to the emitting fluid elements, the 
phase factor becomes 

gifo - ! — cJ2 ) ■ (pi — P2 ) gi^tT (cosh Qi — cosh(y2 —771 )+cosh #2 — cosh(j/i —r/2 )) gZmtr((l — cosh a) (cosh (?i+cosh #2)— sinh a(sinh 8\ — sinh #2)) 



_ iam t T(9 2 ~e 1 -a-a(el+8l)/4)+0(a 3 ,8l 2 ) 



(3.2) 



In anticipation of the saddle-point integration over the flow rapidities to be performed below, I expanded the exponent 
up to second order in 9\ and 2 . Also, I shall be interested mainly in the behavior of the correlation function at small 
relative rapidity and so expanded up to second order in a. 
The same manipulations lead to 



do 1 • K = 



TTTlt 



[(1 + cosh a) cosh 9\ + sinh a sinh 9{\ d6\ d 2 r\ 



of a 
1 + T + 2" 



of a 

1 + T-2" 



dQ 1 d 2 r 1 , 



/ [u\ • e~ (cosh #i+cosh(y2— Vl)) g— [(1+cosh a) cosh 0\ +sinh a sinh 8{\ 

= e --3T-[( 2 + Q2 /2)(l+e 2 1 / 2 )+ae 1 ]+0(a 3 ,6 3 ) 

The expressions for d(T2 ■ A and f e q(u 2 • K) can be obtained by substituting 6\ — > #2 and a - 
the rapidity #1 gives 



(3.3) 
(3.4) 

-a. Integrating over 



ia 2 m t T/ 2 I ^ e -iam t r(e 1 +aeJ/4) e -^(Qe 1 + (l+Q 2 /4)e?) 



a a 
1 + T + 2 " 



a a 
1 + --2" 



'27tT 



1 - ^- (1 + 2?>T) ) e -^^'- J/ 



2 + — ) cos(arT) . 



(3.5) 



On the right-hand side, I performed the following approximations. I expanded the arguments of the exponentials 
assuming a <C 1 but am t /T and oitT of order one. I kept only terms up to 0(1) in this approximation scheme. 
Nevertheless, I dropped a term (i/4)aSm t /T)(TT) 3 = 0(1) since I will mainly be interested in the curvature of the 
correlation function at a = 0. Since (3.5) is symmetric in a — > —a, the integral over 9 2 gives the same result. 
Next, I need Sf(ui ■ K) in the limit 6\ — ► 0, up to second order in a. The result is 



(S/(ui • K) = T c m t 



r 47Y 



(3.6) 



The result for Sf(u 2 ■ K) is the same, except that the derivative of 1/T with respect to rapidity is evaluated at y 2l 
respectively, and a — > —a. Then, 



R(m u a) = J d 2 nd 2 r 2 T 2 2nTm t e - 2mt/T cos 2 {aTT) e^' ^ 5 



X < 1 + T c m t 



1 + tH + £ 3 » 



nil 



T c m t 



T 4Tt ) 

For any function f(y), to leading order in a — y 2 — y\ 

f'(v2) = r(yi) + af"( yi ) . 



1 + T) dT -^ 



I _ " 2 ] 

T ~ 47V J 



(3.7) 



(3.8) 
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For f(y) = 1/T(y) then, 



■ T2 d yi T 



^ logT) 2 - logT - 



(3.9) 



On the right-hand-side, T is evaluated at y±. Integrating over the transverse area, and dividing by the product of the 
single-particle distributions leads to 



C2 [m t , a) — 1 = cos 2 (arT) e 



-rata t T 



I T c m t W 

x < 1 

1 t T 2 



91ogT\ _ 1 l d 2 T 
dlogr) ~ T \dyj 




(3.10) 



Here, (■) refer to averages over the transverse plane as well as over events. For consistency, only terms to first order 
in t c have been kept. 

We can define a longitudinal homogeneity length as the curvature of the correlation function at a = 0: 



1 ld 2 C 2 (a) 



m 2 2 da 2 



t T t T 



rn, 



t 2Tm t { 7 T\dy 2 



'd\ogT 




(3.11) 



The 1 /m 2 prefactor arises because one should in fact take the derivative with respect to the longitudinal momentum 
diff erenc e; for small relative rapidity a, and in the logitudinally comoving frame, P\\/mt = tanha ~ a. The first term 
in (3.11) is sublcading in T/m t <C 1 and will be dropped. Then, 




T 

mi 



T c 1 - 7 ~ Ct + 2A 2 



2r 2 T 2 



(3.12) 



with 



A 2 T 



'dlogT 
, dyi 



T 



d 2 T 
dy\ 



(3.13) 



Deviations from equilibrium manifest themselves in the second term under the square root. They do not distort the 
scaling of J?m with 1/^/nit predicted from a local equilibrium distribution Q (regarding scaling of the out and side 
radii see pp|). No te, however, the "wrong sign" of the correction relative to the leading term, to which we shall return 
below in eq. ( |3.14| ). 

A 2 ^ is the mean-square fluctuation of the temperature^], divided by the typical rapidity scale Ay on which it occurs; 
it measures fluctuations about the average homogeneous particle source^ Only fluctuations on scales Ay smaller than 
the relative rapidity a at which the correlation function 02(a) is being probed^ matter. Fluctuations on larger scales 
do not matter as in that case both particles are emitted from the same rapidity "element" and at the same T, so 
AT = 0. In any case, such rapidity fluctuations of T can only give a significant correction to R\\ if Ay <C 1 and if 
AT/T is of order unity. In other words, after the transition, there should be regions of longitudinal extension ~ 1 fm 
(which at time r ~ 10 fm corresponds to a rapidity interval of 0.1) where T is about equal the ordinary hydrodynamic 
freeze-out temperature for a homogeneous fluid, separated by cold regions of similar size with T ~ (see also p2|). 
Nevertheless, temperature fluctuations of order AT/T ~ 1 appear unrealistically large. Moreover, experimentally 



4 The average temperature field was of course assumed to be rapidity independent, (dT/dy) = 0. 

5 These are analogous to the temperature fluctuations of the cosmic microwave background measured by COBE j2lj, of order 
AT/T ~ 1(T 5 . 

6 The single-inclusive distribution (2.9) is not sensitive to At. That is, of course, due to the fact that it measures the 
temperature at only one rapidity but not correlations of T between two rapidities 1/2 and yi . One might be tempted to measure 
the single-particle distribution in very narrow bins in rapidity, and in individual events, such as to trac e an y local disturbance 
of a boost invariant temperature distribution. This is not possible, however, because the saddle-point (2.4) of the distribution 
function has a width y // 2irT/mt in rapidity, which can not be less than ~ 1 for values of m t where the application of the 
collective hydrodynamic theory makes sense. Thus, temperature fluctuations on scales Ay <C 1 get washed out. 
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one usually measures correlations of identical pions, which are the most abundant hadron species. A significant part 
of the pion yield is commonly believed to originate from decays of resonances (p, 77, A, ...) after freeze-out. Those 
resonance decays smear out any temperature fluctuations on scales Ay -C 1, even if those were present before the 



resonance decays. Therefore, realistically it seems that the correction proportional to A T in eq. (3.12), which would 



tend to enhance R11 , is rather small. 

There is also an increase (1 — 7 > 0) of Ru by an amount depending on the temperature gradient in time direction. 
If that gradient 7 is small, which corresponds to a rather viscous fluid, the increase is largest. If the temperature drops 
steeply, 7 is larger and so the correction is smaller. In any case, by considering either isentropic or isoergic expansion, 
i.e. conservation of either the entropy or the energy in a comoving volume element, one obtains 7 ~ 1/4 — 1/3, and 
so 1-7 = 0(1). 

The only term leading to a reduction of J?ii is that proportional to fr, the curvature of the temperature field in 
rapidity. This is also obvious intuitively. Note that this curvature term arises from the first-order correction to local 
equilibrium (it is proportional to r c ). A nonvanishing curvature of T does not contradict the initial assumption that 
(dT/dy) = at midrapidity. Moreover, one could still argue that T is rapidity independent only on scales probed by 
the single-inclusive distribution, say Ay ~ 1, but does exhibit some local curvature on scales Ay <C 1. Even so, it 
doesn't seem likely that the average local curvature, averaged over the transverse plane and events, is big. 

Therefore, in summary, all of the corrections from eq. ( |3.12] ) proportional to r c should be rather small for the 
case of high-energy nuclear collisions, in particular since they are down by a factor 1/(tT) 2 relative to perfect local 
equilibrium: for nuclear collisions, typically r is on the order of 10 fm and T on the order of 100 MeV, so tT ~ 5. 

This brings us to the main point. From kinetic solutions for the evolution of hadrons produced from a QGP fluid it 
was observed that the two-particle correlation functions are sensitive to a "cloud" of late soft hadronic interactions, 
which only minorly affect the single-inclusive distribution (l^J^l. We can understand qualitatively why this is so. 
Above, I argued that deviations from local equilibrium grow as the temperature T decreases and eventually "freeze" 
the single-inclusive distributi on (|2.9| ), despite the smoothly ongoing cooling down to T = 0. 

The correlation function ( 3.10| ) also exhibits a correction term proportional to t c /t. Note, however, that the 



correction vanishes at small relative rapidity (or longitudinal momentu m) o f the pions, a — > ! The characteristic 
relative rapidity over which C2 falls off is set by the exponential in eq. (3.1C), a 2 ~ l/(mtT 2 T). Thus, the first-order 



correction in the curly brackets in (3.10) is of order r c /(r 3 T ), and can only become large when t c /t ~ r T . This is a 



large number for the case of nuclear collisions, and it grows with time. This means that if we are dealing with a nearly 



perfect hadron fluid (t c /t -C 1) with only linear growth of t c in time, as in eq. (2.11), then non-ideal corrections to Ru 
actually decrease with time. This is despite the fact that they are sufficient to actually "freeze" the single-inclusive 
distribution ! 

For to "freeze", i.e. dR^/dr = 0, the relaxation time must evolve as 

r c = (l-7)--2^T 2 r 2 , (3.14) 
r 1 — 7 

where I dropped £t and At for simplicity. Thus, close to equilibrium, i.e. for small t c /t, R\\ can not "freeze" since 
the relaxation time would have to decrease with time, which seems unphysicalQ. Eq. (3.14) shows that there is no 
physical solution (r c > for r > ro) which leads to "freezing" of Ru in a fluid ne ar eq uilibrium, contrary to N(p). 
That is because for f c to stay positive, the first term on the right-hand-side of eq. ( |3.14| ) must overwhelm the second 
one. Ho wever , the first term by itself asks for a rather slow, nearly linear increase r c ~ r 1-7 . In turn, the negative 



term in (3.14) grows almost quadratically with r, and so will eventually take over. The only way out is to stabilize 
i?H through corrections of higher order in t c , which are not taken into account in ( |3 .11 ). This motivates why the HBT 



radii arc much more sensitive to the dissipative evolution of the hadron "fluid" than the single-inclusive distribution 
or quantities derived from it (flow |Q). 

The "HBT puzzle" can therefore be formulated as follows: if hadronization leads to a hot and dense hadron fluid, 
t c /t is small initially^. It will take a long time until the two-particle correlation function "freezes out" because r c 
first has to evolve slowly into the regime where corrections of higher order in r c become important. Consequently, Ru 
will be large. In fact, assuming that the initial temperature of the hadrons is equal to the confinement temperature 
T c ~ 180 ± 20 MeV already leads to significant deviations of the computed HBT radii f?],|2{| to the data {§. The 



7 At least, if one has in mind typical kinetic estimates of relaxation times, see for example |23| . Near a critical point estimates 
of r c might be altered. 

8 Here, "initial" refers to the initial condition for the hadrons, which is set by the hadronization. 
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theoretical calculations do not fail very badly, but in a clear systematic fashion. While the computed ratio of the out 
and side radii, R /R s , increases with the pair transverse mass m ( , the data shows it rather flat. Also, the experimental 
R\\(mt) does seem to roughly follow the 1/y/mt behavior but is systematically lower than the theory at all m t , even 
for the most "optimistic" assumption T c = 160 MeV. 



IV. HOW IT MIGHT WORK... 



If a slow growth of the relaxation time prevents Ru from freezing quickly, how could one then obtain a small 
homogeneity length ? 

It clearly helps to start the final-state hadron kinetic evolution at low temperature, and with big deviations from 
equilibrium. The closer the inital condition for the hadronic evolution is to the "freeze-out" condition, the better. 
The results of Q already showed that the picture improves if the initial temperature of the hadron fluid is lowered, 
but the use of perfect-fluid dynamics as a hadronization scheme implies that on the hadronization surface r c /r was 
small in those studies. As discussed above, the QGP shouldn't hadronize into an ideal fluid of hadrons, but into a 
very viscous one. 

Employing ideal hydrodynamics to model the phase transition, the fluid of hadrons starts out essentially at T c . 
With T c > 160 MeV, one can not describe the HBT data [f| very well |@,^5|: the hadronic fluid is too dense and r c too 
small. One should probably go lower in temperature, which implies a non-equilibrium phase transition. One way to 
obtain it is from parton cascade models [T^Jl^] where one can tune the parton^hadron transition such as to yield a 
cool hadron fluid with initial temperature To below T c and with large velocity and temperature gradients. (Increasing 
the parton-parton cross section relative to lowest-order pQCD estimates presumably does the job by delaying the 
transition into hadrons.) This might be the reason why ref. [|12[ apparently is able to fit the radii at RHIC (in fact, Ru 
even comes out too small). In field-theory language, in turn, one would have to consider spinodal decomposition |2q ] 
if the transition is first-order. If it is a cross over or weak first-order transition, the decay of a condensate which 
saturates the free energy of the deconfined state at T c [|tJ could be a successful approach: while the dcconfmement 
temperature sets the scale for the effective potential, the reheating temperature of the hadrons is determined by the 
decay process of the condensate and by the overall expansion, and so is quite likely less than T c f]. 

In summary, the pion correlation function requires much smaller relaxation rate than the single-particle distribution 
in order to "freeze" . Therefore, HBT interferometry of the QGP "ashes" |2| might reveal some of the bulk properties 
of the QGP hadronization and its non-equilibrium nature. The first-order theory applied here is not capable of 
quantitative predictions but it clearly points at the fact that "freeze-out" of two-particle correlations at small relative 
momentum is tightly related to the dynamics of the confinement transition. 
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